Extent of sensitivity of single photon production to parton distribution functions 
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We have studied the production of single isolated prompt photons in high energy proton-proton 
collisions at the RHIC (\/s= 200 GeV) and the LHC (\/s= 7 TeV) energies within the frame work 
of perturbative QCD upto next-to leading order of strong coupling (as)- We have used five different 
parameterizations of parton distribution function (PDF) starting from the old CTEQ4M to the new 
CTIO distributions and compared our results with the recent single prompt photon data from the 
PHENIX and the CMS collaborations. The prompt photon cross-section is found to be described 
equally well by all the PDF's within the experimental errors at the RHIC and the LHC energies. 
The deviation in the prompt photon yield for different PDF sets is within ± 20% when compared 
to CTEQ4M. The diphoton measurement could be a potential candidate to constrain the gluon 
distribution inside the proton. 
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The production of single prompt photons of large transverse momenta in high energy hadronic collisions constitutes 
unique signal of the interactions between quarks and gluons at short distances [1]. The single prompt photon yield is 
■ expected to be sensitive to parton distribution function (PDF) in general and to gluon distribution in particular of 
the colliding hadron [2]-[8]. It is also considered an essential ingredient to quantify the nuclear modification of direct 
' photon production in the relativistic nucleus-nucleus collisions [9]. The word "prompt" is used to identify the class 
O I of photons which do not come from the decay of large momenta hadrons (e.g., tt^, t], etc.). The basic mechanisms of 
' prompt photon production at the Leading-order (LO) of strong coupling as are (i) quark-gluon Compton scattering 
(qg~>' q7); (ii) quark-anti-quark annihilation (qq— ^ g7) and (iii) bremsstrahlung radiation from the final state parton 
(— I ' (q(g)— >■ q(g)+7) [1]. Because of the the point-like coupling, the photons produced through the processes (i) and (ii) 
'— are called "direct" photons. The direct photons are often calculated using perturbative quantum chromodyanmics 
^ . (pQCD) which successfully explains the production of particles with energy E ^ Aqc_d (where Aqcd ^ 200 MeV). 
^ The photons originated from the process (iii) are called "fragmentation" photons. The evaluation of this contribution 
requires non-perturbative parton to photon fragmentation function as an input along with pQCD. However, at this 
point we must state that the distinction is renormalization scheme dependent and has no physical implication beyond 
\Q LO. A complete Next-to Leading Order (NLO) theory is needed to explain the single and double prompt photon data 
in hadronic collisions [10]. 

^ The momentum fraction probed by the prompt photons of transverse momentum for proton-proton (p+p) 
. collisions at mid-rapidity (y= 0) is given by; x w 2pi^ / y/s where -y/s is the center of mass energy. The photon 
CO ' production is more sensitive to the gluon distribution at low pi^ (i-C, small x) and to the valance quark distribution 
, towards higher pi^ (i.e., large x) for a fixed -i/s. Thus one can hope to determine the gluon distribution g(a;,Q^) 
J> ' unambiguously from the very accurate prompt photon data. 
• '~j ] In this brief report, we want to verify the above assumptions. For this purpose we have considered the isolated 
rN ■ prompt photon production in p+p collisions at the Relativistic Heavy Ion Collider (RHIC) (•\/s= 200 GeV) and 
^ ] the Large Hadron Collider (LHC) iy/s= 7 TeV) energies. The isolated prompt photon cross section is found more 
" " " advantageous than the inclusive prompt photon cross section [11, 12]. The collider experiments always perform isolated 
photon measurement to identify the prompt photon signal from the large background of secondary decay photons. 
In addition, the isolation criterion makes the total photon cross section much less sensitive to the non-perturbative 
fragmentation process enabling us to study the small x behaviour of gluon distribution more precisely. A photon is 
said to be isolated if the accompanying hadronic transverse energy (Eip'*) inside the isolation cone is less than some 
specified value (E!^°), 



E^ad < E^°; inside R, (1) 



where we considered a cone of radius R {=y/ {y — y-f)^ + (0 — ^7)^) in the rapidity (y) and the azimuthal angle (0) 
space about the photon direction. It has been shown in Rcf. [12] that the QCD factorization theorem is valid in all 
orders of perturbation theory with the above general isolation criterion (Eq. 1). Hence the isolated prompt photon 
cross section is a pertubatively well defined quantity like the inclusive case. In experiment, sometimes the isolation 
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FIG. 1: (Color online) Scale dependence of single isolated prompt photon production in p+p collisions at the RHIC (200 GeV) 
and LHC (7 TeV) center of mass energies, calculated with CTIO PDF and BFG-I fragmentation function. The experimental 
data are adopted from the PHENIX [5] and the CMS [14] collaborations. 

energy is expressed in terms of a dimensionless parameter e — EJ^°/pi^. We write the single isolated prompt photon 
production cross section as [11, 12]: 

Ej-^ (2/7'-P7'R') = X! / d.Xif,/A{xi,n}) dX2fj/B{x2,f4) 

/dz 
—Dcl^{z, ^?p) a]:'^%^x{Pc,xi,X2] Zc,R\ lif , iiR, iiF)Q{z - Zc), (2) 

where the /'s are called the parton distribution functions (PDF's) which give the probability for the incoming 
partons i and j carrying momentum fractions xi and X2 inside the hadrons A and B respectively. Dc/^{z, is giving 
the parton to photon fragmentation probability defined ai z — —. When a photon is emitted from "direct" hard 
interaction (i.e. c = 7), the photon fragmentation function reduces to 5{\ — z). Following the isolation criterion (1); 
we find that the fraction z should be > 1/(1 + e) = Zc [12]. Thus for small values of e, the fragmentation process is 
less contributory to the total photon cross section compared to the direct process. This is contrary to the inclusive 
prompt photon production. The fragmentation function Df./-y{z, fi^) appears in Eq. 2 is the same as in case of 
inclusive photons. is the isolated partonic interaction cross section, contains all process upto 0{a^) (for = 7) 
and upto 0{al) (for c 7^ 7) of the reaction i + j — >■ c + X. The dependence of the isolation parameters (R, e) 
are fully encoded into the isolated cross section (0-*^°). fif,iiji,^p are the three different scales associated with the 
factorization, renormalization and fragmentation respectively. 

With this theoretical introduction, now we show the results of isolated prompt photon production for p+p collisions 
at the RHIC and LHC center of mass energies ( Fig.l). Our calculation is based on the Monte Carlo programme 
JETPHOX [13] which evaluates direct and fragmentation photon cross-section separately at NLO accuracy. We sum 
the two contributions to get the total physical cross section at NLO. JETPHOX has been found to explain the isolated 
prompt photon data at the Tevatron energy [3] and at the LHC energy [2, 14]. This programme is very useful in the 
sense that one can define variety of experimental requirements (kinematic, isolation) at the partonic level. Numerous 
sets of PDF are made available in JETPHOX through the LHAPDF library [15]. We have considered isolation cone 
of radius R= 0.4 and constant isolation energy cut in these studies. The parameters of calculations at the RHIC 
are; p^,= 4-20 GeV, \y^\ < 0.5, E!f°= 2 GeV and at the LHC are; p:^= 10-200 GeV, \y^\ < 1.5, Ef°= 5 GeV. We 
have used the BFG-I parton to photon fragmentation function by Bourhis et al. [16] which includes correction in the 
fragmentation function beyond leading logerthemic approximation. The three scales fj,f, fiR, fxp are set equal to a 
common scale /i, to reduce theoretical uncertainties in the calculation. The scale /i further defined as cpi^ with c= 0.5, 
1.0, 2.0. The scale dependence of the isolated prompt photon production, calculated with CTIO PDF [17] is shown 
in Fig.l. The results have been plotted against the prompt photon measurement by the PHENIX collaboration at 
mid-rapidity [5] and by the CMS collaboration in the rapidity range 0.9 < ij <1.44 [14]. Our results quite agree with 
the theoretical curve mentioned in Refs. [5, 14]. The scale /i = l.Opi^ is found to be a better choice though the relative 
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FIG. 2: (Color online) (Left panel) Comparison of differential cross section of isolated prompt photons calculated with CTEQ4M, 
CTEQ6M, CTIO, MRST2001, NNPDFl.O PDF, BFG-I fragmentation function and scale n = l.Op^, for p+p collisions at the 
RHIC energy. (Right panel) The same for p+p collisions at the LHC energy. 
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FIG. 3: (Color Online) The "direct" and "fragmentation" contributions to total isolated prompt photon production in p+p colli- 
sions are displayed for the five different PDF's at the RHIC (left) and the LHC energies (right). The "fragmentation" component 
is scaled by 0.1. 



uncertainty bet'ween the scales is found to be very small. 

Next we have considered four more different parameterizations of PDF namely; CTEQ4M [18], CTEQ6M [19], 
MRST2001 [20] and NNPDFl.O [21]. Our motivation is to check whether the single prompt photon yield is able to 
distinguish between these PDF sets. We do hope that the recent prompt photon measurements at the RHIC and LHC 
experiments can constrain the gluon density in the proton. The invariant prompt photon yield calculated with the five 
PDF's at the RHIC and the LHC energies is displayed in Fig. 2. The PHENIX collaboration at the RHIC measures 
the inclusive prompt photon spectra 'iJ^d'^cr/d'^p^ 'whereas the CMS collaboration at the LHC measures the isolated 
photon spectra '' da / dp7pdy^\ We have normalized the JETPHOX output accordingly to compare with the experimental 
measurements. We see that starting from very old parameterization like; CTEQ4M to a recent parameterization like; 
CTIO gives reasonably good description of the data within the statistical errors. A little wiggling observed at different 
theoretical curves for pi^ > 12 GeV at the RHIC and for pi^ > 120 GeV at the LHC can be attributed to the Monte 
Carlo technique of integration. 

The direct and fragmentation contributions to the isolated photon cross section (da / dp^dyj) at NLO are shown 
individually in Fig. 3. The NLO corrections to the direct contribution is due to additional real and virtual gluon 
emission at the LO processes qg— )• and qq— g7 [22]. The NLO correction to the fragmentation is coming from 
evaluation of generic 0{al) subprocesses; i + j — c + d + e where the parton c fragments into a photon [23]. All these 
higher oder corrections are included in JETPHOX according to the MS renormalization and factorization scheme. It is 
known indeed that the total cross section is only the physical quantity at NLO because a direct 0{a1) process like; 
qg— >■ qg7 contributes the same as LO fragmentation process when the fragmentation scale fip is close to 1 GeV [12]. 

We expect qg (or qg) initiated "direct" photons could be a good probe to track down the varying gluon distribution 
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FIG. 4: (Color Online) The PDF uncertainty in the isolated prompt photon production in p+p collisions related to the choice 
of PDF is plotted with the transverse momentum of photon (p^) at the RHIC (left) and LHC (right) energy. 
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FIG. 5: (Color Online) The gluon momentum density a::g(a::,Q'^) vs the momentum fraction x for the five PDF's considered in 
the present study for a typical value of momentum transfer square 80 GeV^. 



g(x,Q^) among the PDF sets. However as shown in Fig. 3, we have not found any dependence of the direct and 
fragmentation contributions to PDF at the RHIC and the LHC energies. The fragmentation contribution shown in 
Fig. 3, is scaled down by a factor of 10 for better visibility. We find the direct component has dominant contribution 
(~ 80-90%) to the total photon cross section at the RHIC and the LHC energies. 

Still now we have plotted our results on logerthemic scale which can easily consume a factor of 2-3 difference 
among the PDF's. To measure the deviation in the isolated prompt photon yield of different PDF sets, we introduce a 
parameter Apdf — (Modern — Old) /Old. The "Old" refers to the isolated prompt photon cross section calculated with 
CTEQ4M. The "Modern" refers to the same calculated with CTEQ6M, CTIO, MRST2001 and NNPDFl.O. Thus, the 
CTEQ4M is chosen as central value and Apdf measures the deviation as we go towards newer structure functions. 
The dependence oi ApDp with at the RHIC and the LHC energies, is displayed in Fig. 4. Apart from the numerical 
fluctuations, we find all the PDF's show a monotonic variation with pi^ e.g., Apdf of CTEQ6M first decreases for 
prpT < 10 GeV at the RHIC and < 40 GeV at the LHC, then increases with pi^. The amount of deviation seen 
is least for MRST2001 and most for NNPDFl.O. It is found that the uncertainty in the isolated photon cross section 
associated with the choice of PDF's is within 20% at the RHIC and LHC energy for all p^. Similar work at the 
Tevatron energy can be found in Ref. [24] where the comparison is made with three PDF's and prediction is given for 
p+p collisions at the top LHC energy {^/s= 14 TeV). 

Finally, we have checked the variation of gluon distribution g(x,Q^) among the five different parameterizations of 
PDF in the kinematic regime, 10~^ ^ x < 10~^, probed by the present RHIC and LHC photon data. In Fig. 5, we 
have plotted the momentum density of gluons a;g(a;,Q^) vs the momentum fraction x for a typical value of Q^= 80 
GeV^. This shows that the gluon distribution is varying indeed among the five PDF's considered in this study. We 
believe that the sensitivity of single prompt photon production to the parton distributions is taken away due to the 
integration over the initial parton momenta xi and X2 in Eq. 2. 
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However, the double prompt photon production cross section in hadronic coUisions, is directly proportional to the 
incoming parton momentum density (or the PDF) [25]. The double prompt photons produced through quark-anti- 
quark annihilation qq— ^ 77 or gluon-gluon fusion gg— >■ 77 completely specify parton momentum distribution. The 
single bremsstrahlung contribution like; qg— >■ q77 to the double photon production is also sensitive to the PDF. Hence 
we hope the double prompt photon measurement at the LHC energy could better probe the gluon distribution g(a;,Q^) 
in the proton in comparison with the single prompt photon measurement. 

In summary, we have analyzed the data of prompt photon production cross section in p+p collisions at the center of 
mass energies 200 GeV (RHIC) and 7 TeV (LHC) with five different parameterizations of PDF. Our calculation relies 
on the available programme JETPHOX which calculates isolated prompt photon production at the NLO accuracy. We 
show that the single isolated prompt photon spectra is almost insensitive to the choice of PDF. The present analysis 
suggests the recent prompt photon measurement at RHIC and LHC can constrain the gluon distribution in the proton 
not more than 20% of certainty. We expect the differential cross section of double prompt photons can determine the 
gluon distribution more precisely and shall be checked in future studies. 

Acknowledgements: SD is thankful to Dincsh Srivastava, Rupa Chatterjee for many stimulating discussions and 
comments on the manuscript. The author sincerely acknowledges the correspondences with Jean-Philippe Guillet 
regarding the set-up of JETPHOX. The author is financially supported by Department of Atomic Energy, India during 
the course of this work. 



[1] J. F. Owens, Rev. Mod. Phys. 59, 465 (1987). 

[2] G. Aad et al. [ATLAS Collaboration], Phys. Rev. D 83, 052005 (2011). 
[3] T. Aaltonen et al. [CDF Collaboration], Phys. Rev. D 80, 111106 (2009); 

V. M. Abazov et al. [DO Collaboration], Phys. Lett. B 639, 151 (2006) [Erratum-ibid. B 658, 285 (2008)]. 
[4] S. S. Adler et al. [PHENIX Collaboration], Phys. Rev. Lett. 98, 012002 (2007). 
[5] A. Adare et al. [PHENIX Collaboration], Phys. Rev. D 86, 072008 (2012). 

[6] P. Aurenche, R. Baier, M. Fontannaz, J. F. Owens and M. Werlen, Phys. Rev. D 39, 3275 (1989). 

[7] H. Baer, J. Ohnemus and J. F. Owens, Phys. Rev. D 42, 61 (1990). 

[8] W. Vogelsang and A. Vogt, Nucl. Phys. B 453, 334 (1995). 

[9] F. Arleo, JHEP 0609, 015 (2006). 
[10] W. Vogelsang and M. R. Whalley, J. Phys. G 23, Al (1997). 
[11] E. L. Berger and J. Qiu, Phys. Rev. D 44, 2002 (1991). 

[12] S. Catani, M. Fontannaz, J. P. Guillet, and E. Pilon, JHEP 0205, 028 (2002). 
[13] http : //lapth . cnr s . f r/PHOX_FAMILY/ j etphox . html; 

P. Aurenche, J. P. Guillet, E. Pilon, M. Werlen, and M. Fontannaz, Phys. Rev. D 73, 094007 (2006). 
[14] S. Chatrchyan et al. [CMS Collaboration], Phys. Rev. D 84, 052011 (2011). 
[15] http : //www.hepf orge . org/downloads/lhapdf . 

[16] L. Bourhis, M. Fontannaz and J. P. Guillet, Eur. Phys. J. C 2, 529 (1998). 

[17] H. L. Lai et al, Phys. Rev. D 82, 074024 (2010). 

[18] H. L. Lai et al, Phys. Rev. D 55, 1280 (1997). 

[19] J. Pumplin et al, JHEP 0207, 012 (2002). 

[20] A. D. Martin et al, Eur. Phys. J. C 23, 73 (2002); 

A. D. Martin et al, Eur. Phys. J. C 28, 455 (2003). 
[21] R. D. Ball et al [NNPDF Collaboration], Nucl. Phys. B 809, 1 (2009) [Erratum-ibid. B 816, 293 (2009)]. 
[22] P. Aurenche, R. Baier, M. Fontannaz and D. Schiff, Nucl. Phys. B 297, 661 (1988). 
[23] F. Aversa, P. Chiappetta, M. Greco and J. P. Guillet, Phys. Lett. B 210, 225 (1988); 

P. Aurenche, P. Chiappetta, M. Fontannaz, J. P. Guillet and E. Pilon, Nucl. Phys. B 399, 34 (1993). 
[24] R. Ichou and D. d'Enterria, Phys. Rev. D 82, 014015 (2010). 
[25] E. L. Berger, E. Braaten and R. D. Field, Nucl. Phys. B 239, 52 (1984). 

R. D. Field, Applications of Perturbative QCD, Addison-Wesley Publishing Company. 



